Multipolar correlations and deformation effect on nuclear transi- 
tion matrix elements of double-/? decay 
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Abstract. - The two neutrino and neutrinoless double beta decay of ^*'^^Zr, ^^'^'^''Mo, ^"''Ru, 
"°Pd, i28,i30rpg isoj^jj isotopes for the 0+ 0+ transition is studied within the PHFB frame- 
work along with an effective two-body interaction consisting of pairing, quadrupole-quadrupole 
and hexadecapole-hexadecapole correlations. It is found that the effect of hexadecapolar correla- 
tions can be assimilated substantially as a renormalization of the quadrupole-quadrupole interac- 
tion. The effect of deformation on nuclear transition matrix elements is investigated by varying 
the strength of quadrupolar correlations in the parent and daughter nuclei independently. The 
variation of the nuclear transition matrix elements as a function of the difference in deformation 
parameters of parent and daughter nuclei reveals that in general, the former tend to be maximum 
for equal deformation and they decrease as the difference in deformation parameters increases, 
exhibiting a very similar trend for the (/3~/3~) and (/3~/3~) transition matrix elements. 



The nuclear /3/3 decay is a rare second order semilep- 
tonic transition between two even Z-even N isobars "^X 
and z±2^ involving strangeness conserving weak charged 
currents. It is expected that the nuclear decay can 
proceed through sixteen experimentally distinguishable 
modes, namely double-electron emission double- 
positron emission {f3~^f3~^), electron-positron conversion 
(e/S"*") and double-electron capture (ee) with the emission 
of two neutrinos, no neutrinos, single Majoron and double 
Majorons. The /3+/3^, e/3"'" and ee processes are energeti- 
cally competing and we shall refer to them as e~^P/3 modes. 
The {(}(3)2v decay conserves the lepton number and is an 
allowed process in the standard model of electroweak uni- 
fication (SM). On the other hand, the {(3(3)^^ decay vio- 
lates the conservation of lepton number and can occur in 
a number of gauge theoretical models, namely GUTs -left- 
right symmetric models and E(6)-, Rp-conserving as well 
as violating SUSY models, in the scenarios of leptoquark 
exchange, existence of heavy sterile neutrino, composite- 
ness and Majoron models. Hence, it is a convenient tool 
to test the physics beyond the SM. 



The occurrence of {l3~l3~) 21, decay was proposed by 
Mayer in 1935 [1] as the possibility to explain the long 
observed half-lives of some even-even nuclei and it was 
finally confirmed by Elliott et al. [2] using a TPC detec- 
tor. Presently, the {P~(3~)2i, decay has been already ob- 
served experimentally in ten nuclei, out of 35 possible can- 
didates [3]. For the {e~^/3/3)2^ modes, experimental limits 
on half-lives have been also given for 24 out of 34 possible 
isotopes [3]. Klapdor-Kleingrothaus and his group have 
reported the observation of (/3~/3~)oi/ decay of ^^Ge with 
T°^^ = 2.23j:^:3i X 10^5 yr [4,5]. The first results [4] were 
debated [6-9]. These results could be confirmed by the 
forthcoming MAJORANA [10] and GERDA [11] experi- 
ments. It would be very exciting if this very first obser- 
vation of {f3~P~)oi, decay by Klapdor-Kleingrothaus and 
his coworkers is verified by the latter two mentioned ex- 
periments. 

It is possible to extract accurate nuclear transition ma- 
trix elements (NTMEs) M21, from the observed half-lives 
of {l3~l3~)2i, decay for the 0+ 0+ transition and a com- 
parison between the theoretically calculated and experi- 
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mentally extracted NTMEs provides a cross-check on the 
reliability of different nuclear models used for the study of 
nuclear /3/3 decay. Using the average half-lives [12] and the 
phase space factors [13], the extracted NTMEs M21, vary 
from 0.0152±0.0008 (0.0238±0.0013) to 0.1222±0.0034 
(0.1909±0.0053) for gA = 1.25 (1.00) corresponding to 
^^°Te and ^°°Mo respectively. In all cases, it is observed 
that the NTMEs are sufficiently quenched. Vogel and 
Zirnbauer were the first to provide an explanation of the 
observed suppression of M21/ in the QRPA employing a 
schematic siirfacc-delta interaction including the particle- 
particle channel, whose strength controls both the con- 
centration of P~ strength in the GTR and the reduction 
of total /?+ strength and its concentration at low exci- 
tation energies [14]. These conclusions were verified by 
Civitarese et al. [15] in the same QRPA model using re- 
alistic interaction derived by G-matrix method. In the 
case of {(}~ (3~)ov decay, the nuclear models usually predict 
half-lives assuming certain value for the neutrino mass, or 
conversely extract various parameters from the observed 
limits on half-lives. Over the past years, many a nuclear 
models have been employed to study the nuclear de- 
cay [16,17]. The experimental as well as theoretical study 
of nuclear (3(3 decay is quite wide in scope and has been 
excellently reviewed over the past decades, which can be 
found in the recent review [18] and references there in. 

The structure of a nucleus is mainly decided by the 
subtle interplay of pairing and quadrupolar correlations 
present in the effective two-body interaction. The former 
is responsible for the sphericity of the nucleus, whereas 
the latter makes the nucleus deformed. In addition to the 
pairing interaction, which plays an important role in all 
even Z-even N (3(3 emitters, it has been also shown that 
the deformation degrees of freedom are crucial in the struc- 
ture of i""Mo and ^^^Nd isotopes [19]. The shell model is 
the best choice to calculate the NTMEs as it attempts to 
solve the nuclear many-body problem as exactly as pos- 
sible. On the other hand, the QRPA and its extensions 
have emerged as the most successful models in correlating 
single-/? GT strengths and half-lives of {(3~P~)ov decay. In 
spite of the spectacular success of the QRPA in the study 
of PP decay, there was a need to include the deformation 
degrees of freedom in its formalism. 

Over the past years, the deformed QRPA model has 
been developed for studying P(3 decay of spherical as well 
as deformed nuclei. The comparison of the experimen- 
tal GT strength distribution B{GT) with the results of 
QRPA calculations was employed as a novel method of 
deducing the deformation of '^^Sr nucleus [21]. The ef- 
fect of deformation on the {(3~ (3~)2u decay for the ground 
state transition "^^Ge ^^Se was studied in the frame- 
work of the deformed QRPA with separable GT residual 
interaction [22] . A deformed QRPA formalism to describe 
simultaneously the energy distributions of the single-/3 
GT strength and the {[3~ (3~)2 decay matrix elements 
of 48Ca, 76Ge, 82Se, ^^Zr, i"OMo, "ecd, 128,130^0, i36xc 
and ^^°Nd isotopes using deformed Woods-Saxon poten- 



tials and deformed Skyrme Hartree-Fock mean fields was 
developed [23]. Rodin and Faessler [24] have studied the 
(3(3 decay of ^^Ge, i°°Mo and ^^^Te isotopes and it has 
been reported that the effect of continuum on the NTMEs 
of {(3~(3~)2^, decay is negligible whereas the NTMEs of 
{(3~ (3~)q^ decay are regularly suppressed. 

Recently, large scale calculations in the "Interacting 
Shell Model" (ISM) have been performed to study the de- 
pendence of NTMEs for (3- (3- decay of ^^G&, ^^Ge, ^^Se, 
iiopd, ii^Cd, i24Sn, i28,i30Te and ^^exe nuclei on the va- 
lence space, effective two-body interactions and deforma- 
tion [25] . In the ISM, the effect of pairing and quadrupolar 
correlations on the NTMEs of {P~(3~)q^, decay has been 
also studied, in which it has been noticed that the NTMEs 
are the largest in the presence of former only [26,27]. It 
has been shown by Menendez et al. [27] that as the differ- 
ence of deformation between parent and grand daughter 
grows, the NTMEs for both {(3' (3')^^ and (/3"/3")o^ de- 
cay decrease rapidly. 

The PHFB model in conjunction with pairing plus 
quadrupole-quadrupole (PQQ) interaction [28] is a conve- 
nient choice to study the role of quadrupolar correlations 
vis-a-vis deformation on the NTMEs for /3/9 decay. How- 
ever, the PHFB model is unable to provide information 
about the structure of the intermediate odd-odd nuclei in 
its present version and hence, on the single /3 decay rates 
and the distribution of GT strength. In spite of this limita- 
tion, the PHFB model in conjunction with PQQ interac- 
tion has been successfully applied to study the {(3~ (3~)2^ 
as well as {e'^ (3(3)^^ decay modes for the 0+ 0+ transi- 
tion where it was possible to describe the lowest excited 
states of the parent and daughter nuclei along with their 
electromagnetic transition strengths, as well as to repro- 
duce their measured decay rates [29-31]. 

An inverse correlation between the GT strength and 
quadrupole moment has been noticed [20]. In the PHFB 
model, the existence of an inverse correlation between the 
quadrupole deformation and the size of NTMEs M2,y has 
been also confirmed [29-31]. In addition, it has been 
observed that the NTMEs are usually large in the ab- 
sence of quadrupolar correlations. With the inclusion of 
the quadrupolar correlations, the NTMEs are almost con- 
stant for small admixture of the QQ interaction and sup- 
pressed substantially in realistic situation. In the mass 
range A = 90 — 150, the suppression factor for NTMEs 
evaluated without and with the quadrupole-quadrupole in- 
teraction ranges between 2 and 6, both for the (/?"/?") 
and {(3~ (3~ )qj, decays. For the double beta decay of ^^"Nd, 
it was also found [35] that the NTMEs for (/3"/3")o^ de- 
cay have a well defined maximum when the deformation 
of parent and daughter nuclei are similar and they are 
suppressed for a difference in the deformation. 

The purpose of the present work is two fold. In the 
first part, we investigate the effect of hexadecapolar cor- 
relations on the previously calculated spectroscopic prop- 
erties, NTMEs M2y for (/3"/3")2„ [29, 31] and M^^ for 
{(3~ (3~)qjj decay [35]. In the second part, we investigate 
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the effect of having different deformations in parent and 
daughter nuclei on NTMEs M2,y and M^v of ifi^ P^)2i, ^iid 
{(3~P~)q^ decay respectively for the 0+ — > 0+ transition. 

The inverse half-life of the {0~ (i~)2v decay for the 0+ — > 
0+ transition is given by [32-34] 



with 



[Ti%(0+^0+)]-i=G2.|M2.| 



(1) 



The phase space factor G2v can be calculated with good 
accuracy [13,33] and the nuclear model dependent NTME 
is written as 



M2u 



= E 

N 

= E 



(0^lkr+||l+)(l+||aT+||0+) 
En - {Ei + Ef)/2 

(0+||ar+||l+)(l+||aT+||0+) 



with 
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Eq + En — Ej 



Ef) = ^Ql30 + me 



(2) 



(3) 



The expressions to calculate in the PHFB model are 
extensively discussed in our earlier works [29,31,35]. In 
what follows, we include only the required extensions in 
the formalism. The HFB wave functions are generated 
using an effective Hamiltonian with PQQHH type of ef- 
fective two-body interaction. Explicitly, the Hamiltonian 
can be written as 

H = Hsp + V{P) + Cgg [V{QQ) + V{HH)] (4) 

where Hgp denotes the single particle Hamiltonian. 

The V{P), V{QQ) and V{HH) represent the pairing, 
quadrupole-quadrupole and hexadecapole-hexadecapole 
part of the effective two-body interaction. The is 
an arbitrary parameter and the final results are obtained 
by setting the Cqq = I- The purpose of introducing is 
to study the role of deformation by varying the strength 
of QQHH interaction. 

The pairing part of the effective two-body interaction 
V{P) is given by 



V{P) 



alala^af} (5) 



a/3 



where a denotes the quantum numbers {nljm) and the 
state a is same as a but with the sign of m reversed. The 
QQ part of the effective interaction V{QQ) is expressed 
as 

V{QQ) = - J2 E("^)''^"I«2mI7>(/3|92-^|5) at 



where 
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The HH part of the effective interaction V{HH) is given 
as 



V{HH) 



Y.Y. 



(8) 



(9) 



The relative magnitudes of the parameters of the HH part 
of the two body interaction are calculated from a relation 
suggested by Bohr and Mottelson [39]. According to them 
the approximate magnitude of these constants for isospin 
T = is given by 



47r 



XA = 



2X + IA (r 



,2A-2\ 



/orA = 1,2,3,4. 



(10) 



and the parameters for the T = 1 case arc approximately 
half of their T = counterparts. The value of Xi for T = 1 
is taken as exactly half of the T = case. Explicitely 



X4 



2/3' 



0.4884x2-4-2/36" 



(11) 



with h = 1.0032^1/6. 

The model space, single particle energies (SPE's), pa- 
rameters of PQQ type of effective two-body interactions 
and the method to fix them are the same as in refer- 
ences [29,31,35]. We fix X2pn as well as Xipn through 
the experimentally available energy spectra for a given 
model space, SPE's, Gp, Gn, X2pp and xapp- The strength 
of proton-neutron [pn) component of the QQ interaction 
X2pn is varied so as to obtain the spectra of considered nu- 
clei, namely ^^^^^z^^ 94,96,98,100^0^ 98,ioo,io4j^^^ io4,iiopd^ 

iiocd, i28,i30Te, I28,i30xe, isoNd and i5°Sm in optimum 
agreement with the experimental results. The theoretical 
spectra has been taken to be the optimum one if the exci- 
tation energy of the 2+ state E2+ is reproduced as closely 
as possible to the experimental value. All these input pa- 
rameters are kept fixed to calculate other spectroscopic 
properties. 

In table 1, we present the values of X2pn for the PQQ 
and PQQHH interactions. It is clear from the fitted val- 
ues of X2pn that once the procedure to fix E2+ is adopted, 
the change in the latter is minimal even in the presence of 
HH correlations. This behavior is further reflected in the 
calculated yrast spectra, reduced B{E2:Q'^ — > 2+) transi- 
tion probabilities, static quadrupole moments (5(2+) and 
gyromagnetic factors 17(2"'"), which are usually in an overall 
agreement with the experimental data [29,31]. The maxi- 
mum change in all the calculated spectroscopic properties 
is about 10% except for the B{E2:Q+ 2+) and (3(2+) of 
^^Zv, which change by 20% and 23% respectively. Further, 
we tabulate the theoretically calculated and experimen- 
tal deformation parameter /32 [36] of parent and daugh- 
ter nuclei involved in /3^/3^ decay in the same table 1 at 
Cqq = 1. The calculated NTMEs M2u for the 0+ 0+ 
transition at (qq = 1 are given in table 2. The change 
in results for P2 are below 10%. The NTMEs M21, also 
change up to 10% except ^^Zr, ^^OTe and ^^^Nd, for which 
the changes are 17%, 26% and 21% respectively. The ra- 
tio for deformation effect D2^ are 2.29(2.23), 3.70(3.38), 
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Table 1: Fitted values of Xzpn, calculated deformation param- 
eters f32 of parent(P) and daughter(D) nuclei participating in 
l3~ P~ decay using PQQ and PQQHH type of two-body in- 
teraction along with the average experimental values, f and % 
denote PQQ and PQQHH interaction respectively. 
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1.86(1.75), 2.33(2.09), 5.47(5.19), 3.14(3.30), 4.26(3.97), 
2.89(3.31) and 5.94(6.46) for 94,96zj., ss.ioojyjQ^ ^O'^Ru, 
"Opd, i28,i30Te and ^'^"m nuclei with PQQ{PQQHH) 
interaction respectively and their variation due to the 
PQQHH type of interaction is up to 15%. 

A suppression of NTMEs for 0/3 decay with respect 

to the spherical case has been reported when the parent 
and daughter nuclei have different deformations [23,27]. 
To investigate this effect, we present the NTMEs for the 
(/3-/?-)2„ as well as (/?-/?-)o^ decay of^^^^^Zi, ^»^^°"Mo, 
lO'^Ru, iiopd isotopes in fig. 1 and i28,i30rpg 150-^^ 



Table 2: Calculated NTMEs Ma^ and Mq^ for {P' P') 

and (/3~/3~)g decay respectively using PQQ [29,31,35] and 
PQQHH interaction along with the extracted M2v from av- 
erage/recommended experimental half-lives for the 0^ 0"*" 
transition [12], except for ^*Zr [37], and "°Pd [38]. The num- 
bers corresponding to (a) and (b) arc calculated for qa = 1.25 
and 1.0 respectively, f and J denote PQQ and PQQHH inter- 
action respectively. 
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isotopes in fig. 2 as a function of the difference in the 
deformation parameter A/32 = P2{parent) — P^idaughter) 
between the parent and daughter nuclei. The NTMEs are 
calculated by keeping the deformation for parent nuclei 
fixed at Qqq = 1 and changing the deformation of daughter 
nuclei by varying Qqq in the range 0.0 - 1.5.^ 

In our earlier works [29,31], it has been shown that the 
NTMEs are usually large for C,qq = 0.0 i.e. when both 
the parent and daughter nuclei are spherical. With the in- 
crease of (qq, the NTMEs remain almost constant and then 
decrease around the physical value (qq = 1.0 establishing 
an inverse correlation between and 02- It can be ob- 
served at the left panels of figs. 1 and 2 that in all cases 
that when the absolute value of the difference in deforma- 
tion is close to zero the largest NTMEs are obtained. In 
addition, the effect of hexadecapolar correlations is almost 
negligible. 

In the mass mechanism neglecting pseudo-scalar and 
weak magnetism terms of the recoil current, the impor- 
tance of which has been reported by Simkovic et al. [40,41] 
and ncx'ds furtlu^r investigation, the inverse half-life of the 

-"-In table 5 of rcf. [.31], the given 132 values of ^^"Xc (for 
iSOxe^iaOxe) at Cgg = 0.9 and 1.5 are 0.159 and 2.155 respec- 
tively, which is a typographical error. The correct values are 0.157 
and 0.216 respectively. 
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Fig. 1: NTMEs of «^'»«Zr, 9«'i°°Mo, i°^Ru, ""Pd isotopes for 
{p~fi~^^ (left hand side) and (right hand side) de- 

cay as a function of the difference in the deformation parameter 
A/32, "x" denotes the value of NTME for calculated A/32 from 
table 1. 



(/3 /3 )qj^ decay due to the exchange of light neutrinos for 
the 0"*" 0+ transition is given by [32-34] 
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Fig. 2: NTMEs of i28.i30Tg izo^^ j^^^jg^ f^j. ^^-^-•^^^ 
and ) decay. Further details are given in fig. 1. 



where the NTME Mq^, is given by 
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The calculation of Mqi, in the PHFB model has been 

discussed in an earlier work [35]. It is worth mention- 
ing that we include the finite size and short range ef- 
fects through a dipole form factor and Jastrow type of 
short range correlation (SRC). The other two prescrip- 
tions for including SRC namely, through the exchange 
of p and u> mesons [42] and unitary correlation operator 
method [43, 44] need to be investigated. In table 2, we 
give the theoretically calculated Mq^ using the HFB wave 
functions in conjunction with PQQ [35] and PQQHH 
interaction. The change in the values of NTMEs Mo,y 
are below 10% except ^"'Pd, "°Te and ^^ONd, for which 
the changes are 17%, 18% and 25% respectively. The ra- 
tio for deformation effect Do,^ are 2.51(2.30), 4.50(4.21), 
1.94(1.87), 2.15(1.92), 3.81(3.85), 2.64(2.87), 4.43(3.85), 
2.96(3.26) and 6.16(7.14) for ^^'^^Zr, 98,100^0^ i04j^y^ 
"Opd, i28,i30rp^ i50N(J jj^cigj ^jth PQQiPQQHH) 
interaction respectively corresponding to a maximum 
change up to 16%. It is worth mentioning that, the 02 
for ^^°Nd is close to the experimental value whereas for 
^^°Sm, it is off. For the PQQ interaction, the experimen- 
tal /32 but for other physical properties can be reproduced 
with (qq = 0.96 and the corresponding Mq,^ is 0.7015. 

In the right panel of figs. 1 and 2, we present the vari- 
ation of jMoi/] with respect to |A/32| for the above men- 
tioned (3~ (3~ emitters. It is noticed that variation in \Moiy\ 
with changing |A/32| is similar to that of {P~P^)2i, decay. 
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Moreover, it is noticed from figs. 1 and 2 that the NTMEs 
remain constant even when one of the; nuclei is splic;rical 
or slightly deformed. With further increase in deforma- 
tion, the NTMEs become the maximum for |A/32| = 
and then decrease with increase in the difference between 
the deformation parameter. 

We conclude that the eff'ect of including the hexade- 
capole interaction channel is minor both on the deforma- 
tion (^2 and on the NTMEs. In particular, it is clear that 
the qualitative dependence of the NTMEs on deformation 
is not affected by including the HH interactions in the 
Hamiltonian, and the changes, while visible, are small. It 
exhibits that the inclusion of QQ interaction is a neces- 
sary ingredient for a reliable calculation of NTMEs, and 
that effect of the hexadecapolar correlations can be mim- 
icked by a slight renormalization in the QQ channel. The 
independent deformations of initial and final nuclei are, 
indeed, crucial parameters to describe realistic NTMEs 
M21, and Mqv of /3(3 decay, which exhibit, in their respec- 
tive scales, a very similar dependence on the difference of 
deformations. 
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